Evolution of order and chaos across a first-order quantum phase transition 
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We study the evolution of the dynamics across a generic first order quantum phase transition in an interacting boson model of nuclei. 
The dynamics inside the phase coexistence region exhibits a very simple pattern. A classical analysis reveals a robustly regular 
dynamics confined to the deformed region and well separated from a chaotic dynamics ascribed to the spherical region. A quantum 
analysis discloses regular bands of states in the deformed region, which persist to energies well above the phase-separating barrier, 
in the face of a complicated environment. The impact of kinetic collective rotational terms on this intricate interplay of order and 
chaos is investigated. 
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Quantum phase transitions (QPTs) are qualitative changes in 
1 the properties of a physical system induced by a variation of 
parameters A in the quantum Hamiltonian H(A) Such 
ground-state transformations have received considerable atten- 
tion in recent years and have found a variety of applications in 
many areas of physics and chemistry |3fl. The competing inter- 
i actions in the Hamiltonian that drive these transitions, can af- 
fect dramatically the nature of the dynamics and, in some cases, 
lead to the emergence of quantum chaos. This effect has been 
observed in quantum optics models of N two-level atoms inter- 
acting with a single-mode radiation field |4J], where the onset 
of chaos is triggered by continuous QPTs. In this case, the un- 
derlying mean-field (Landau) potential V{A) has a single mini- 
mum which evolves continuously into another minimum. The 
situation is more complex for discontinuous (first-order) QPTs. 
Here V(A) develops multiple minima that coexist in a range of 
A values and cross at the critical point, A = A c . Understanding 
the nature of the underlying dynamics in such circumstances is 
a primary goal of the present Letter. 

The interest in first-order QPTs stems from their key role 
in phase-coexistence phenomena at zero temperature. In con- 
densed matter physics, it has been recently recognized that, for 
clean samples, the nature of the QPT becomes discontinuous 
as the critical-point is approached. Examples are offered by 
the metal-insulator Mott transition 0], itinerant magnets |0], 
heavy-fermion superconductors J3l, and quantum Hall bilay- 
ers (8D . First-order QPTs are relevant to shape-coexistence in 
mesoscopic systems, such as atomic nuclei J3l, and to optimiza- 
tion problems in quantum computing fioll . 

Hamiltonians describing first-order QPTs are often non- 
integrable, hence their dynamics is mixed. They form a sub- 



class among the family of generic Hamiltonians with a mixed 
phase space, in which regular and chaotic motion coexist. In the 
present Letter, we wish to illuminate, in a transparent manner, 
those aspects of this mixed dynamics which reflect the first- 
order transition. For that purpose, we employ an interacting 
boson model which describes such QPTs between spherical and 
axially-deformed nuclei. Our main results are that, (i) in spite 
of the abrupt structural changes taking place across a first-order 
QPT, the dynamics in the coexistence region exhibits a very 
simple pattern. A robustly regular dynamics is confined to the 
deformed region, and is well separated from the chaotic dynam- 
ics ascribed to the spherical region, (ii) The deviations from this 
marked separation is largely due to kinetic collective rotational 
terms in the Hamiltonian. This simple pattern of mixed dynam- 
ics was initially observed at the critical point, A = A C ll 111 . Here 
we show it to be a hallmark of the whole coexistence region. 
Simply divided phase spaces were encountered in billiard sys- 
tems Il2lll3ll . which are generated by the free motion of a point 
particle inside a closed domain whose geometry governs the 
amount of chaoticity. Here, in contrast, we consider many -body 
interacting systems undergoing QPTs, where the onset of chaos 
is governed by a change of coupling constants in the Hamilto- 
nian. 
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The interacting boson model (IBM) [ 14] describes 
quadrupole collective states in nuclei in terms of a system of N 
monopole (s) and quadrupole (d) bosons, representing valence 
nucleon pairs. The Hamiltonian conserves the total boson num- 
ber N and angular momentum L. Its geometric visualization 
is obtained by a potential surface, V(J3, y) - <J3, 7; N\H\/3, y; N), 
defined by the expectation value of the Hamiltonian in the in- 
trinsic condensate state \/3, y; N) = (N\y 1/2 [r f c (f3, y)] N \0 ), where 
rl(J3,y) = ±[fJcosydl+f3sm 7 ±(4+dl 2 )+ V^eV] Q. 
Here (J3, y) are quadrupole shape parameters whose values 
(y0eq,7eq) at the global minimum of V(J3, y) define the equi- 
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librium shape for a given Hamiltonian. QPTs between such 
stable shapes have been studied extensively in the IBM frame- 
work JT^, [13, EM] and are manifested empirically in nuclei 
Their nature is dictated by the topology of V(J3, y), which 
serves as a Landau potential. For that reason, in studying such 
QPTs, it is convenient to resolve the Hamiltonian into two 
parts d], 
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The intrinsic part (//j n t) determines the potential surface V(J3, y), 
while the collective part (H co \) is composed of kinetic terms 
which do not affect the shape of V(J3, y). 

Focusing on first-order QPTs between stable spherical (y6 eq = 
0) and prolate-deformed (J3 eq > 0, y eq = 0) shapes, the intrinsic 
Hamiltonian reads 



Figure 1 : Behavior of the order parameter, /3 cq , as a function of the control pa- 
rameters (p, f) of the Hamiltonian {2). Here p" , (p c , f c ), f*, are the spinodal, 
critical and anti-spinodal points, respectively, given by p* = A* [-(r 2 - 4r+ 1) + 

(r^jfTlFlffl] 1 ' 2 , (Pr =/?(,', & = 0), r* = (l +PlT\ withr=/; 2 . 
The deformation at the global (local) minimum of the Landau potential {3} is 
marked by solid (dashed) lines. /3 cq =0 |y3 c q = V2/3o(1+/Sq) '' 2 ] on the spherical 
[deformed] side, with values shown correspond to /?o = 1.35. 



ff m (p)/h 2 = 2{l-p 2 l3l)h d {h d -l)+l3lRl-R2, (2a) 
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the following classical potential 



where hj - 2^ d^d^ is the if-boson number operator, R 2il (p) = 
^Hs% +p Vf&iW)® , PjOSb)^ ■ d^ -f (s f ) 2 and p\jfi^) = 
yftfosUl + V7(rftrft)(2). Here R ^ = { _iy Rl ^ p 2ii = 
(—lyP^-ii and the dot implies a scalar product. Scaling by 
hi = h 2 IN{N - 1) is used throughout, to facilitate the compar- 
ison with the classical limit. The control parameters that drive 
the QPT are p and £ with < p < /J" 1 and f > 0, while p 
is a constant. For the indicated ranges, the intrinsic Hamilto- 
nians in the spherical [H' nt (p)] and deformed [HEXif)] phases 
have the intrinsic states \p eq ,y eq ;N) with, respectively, p eq = 
and [p eq >0, y eq = 0], as zero-energy ground states. For large N, 
the normal modes of H.ip) [H.JJ^)] involve quadrupole [both 
P and y] vibrations about the spherical [deformed] global min- 
imum, with frequency e = 2h 2 Npl [ep = IhiNp 2 ^!^ + 1), e y = 
18^Vj3q(1 +/3q)~ 1 ]. The two Hamiltonians coincide at the crit- 
ical point p c =0^ 1 and £. = 0: H\ Bt (p c ) = ff^ t (&), being equal to 
the Hamiltonian studied in 111 111 . 

The classical limit of the IBM is obtained through the use 
of coherent states and taking N — > 00, with l/N playing 
the role of h H2CH1 . Number conservation ensures that phase 
space is 10-dimensional and can be phrased in terms of two 
shape (deformation) variables, three orientation (Euler) an- 
gles and their conjugate momenta. The classical Hamiltonian 
obtained involves complicated expressions (including square 
roots) of these variables. Setting all momenta to zero, yields 
the classical potential which is identical to V(J3, y) mentioned 
above. Chaotic properties of the IBM have been studied exten- 
sively 112 in . albeit, with a simplified Hamiltonian, giving rise 
to extremely low barrier and narrow coexistence region. The 
recent identification of IBM Hamiltonians without such restric- 



tions [22] enables, for the first time, a comprehensive analysis 
across a generic first-order QPT 

For the Hamiltonian of Eq. ©, the above procedure yields 
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where F = cos3y. The variables p e [0, y2], y e [0, 2n) 
can be interpreted as polar coordinates in an abstract plane 
parametrized by Cartesian coordinates x = /Jcosy and y = 
jfisiny. The potential V'(p) [V n (^)] has a global spherical [de- 
formed] minimum with, respectively, p eq = [p eq > 0, y eq = 0] . 
At the spinodal point (p*), V'(p) develops an additional local 
deformed minimum, and the two minima become degenerate 
at the critical point p c (or £ c ). The spherical minimum turns 
local in V 7/ (£) for £ > % c and disappears at the anti-spinodal 
point (£**). The order parameter p eq , shown in Fig. [1] is a 
double-valued function in the coexistence region (in-between 
p* and ^**) and a step-function outside it. The potentials 
V(p, y = 0) = V(x,y = 0) for several values of i;,p, are shown 
at the bottom rows of Figs. [2][3] [panels (a)-(e)]. The height of 
the barrier at the critical point is V b = h 2 [l - (1 + /3q) 1/2 ] 2 /2. 
Henceforth, we set Pq — 1 .35 which is a typical value within the 
acceptable range, 1 < /Jo < 1.41, for deformed nuclei. In this 
case, Vbjhi = 0.231 (compared to Vb/fi2 = 0.0018 in previous 
works [21]). 

The classical dynamics of L = vibrations, governed 
by Hi nt <GK can be depicted conveniently via Poincare sec- 
tions [23]. These are shown for selected energies below the 
domain boundary En m = V(p = v2, y), and control parameters, 
p < p c in Fig. [2] and ^ > g c in Fig. [3] For p = 0, the system 
is integrable, with V ! (p = 0) cc p 1 ^ 1 + i(l -/^jS 4 . The sec- 
tions for p = 0.03 in Fig. Q~l show the phase space portrait typ- 
ical of an anharmonic (quartic) oscillator (AO) with two major 
regular islands, weakly perturbed by the small p cos 3y term. 
For small p, V r (p) ~ 1 - p "V^S 3 cos 3y. The derived phase- 
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Figure 2: (Color online). Poincare sections (upper five rows) depicting the 
classical dynamics of H- m Sp) )2at with ti2 = 1, /So = 1.35, for several values 
of p < p c . The bottom row displays the corresponding classical potentials 
V {I) (p) {3a). The five energies, below £i; m = 2/i2, at which the sections were 
calculated consecutively, are indicated by horizontal dashed lines. The Peres 
lattices jx,-, E;}, portraying the quantum dynamics for eigenstates \i) of S! (p) 
with L = and iV = 80, are overlayed on the classical potentials V(x, y = 0). 



Figure 3: (Color online). Same as in Fig.|2]but for the intrinsic Hamiltonian 
6" t (g) <2bJ and classical potential V (//, (f) @b), with f > £ c and E lim = h 2 (2+£). 
Notice in the Peres lattices | at the bottom row, the sequences of regular 

states in the deformed region x ss 1 [also observed in Figs. 2(d)-2(e)]. The 
lowest sequence consists of (TV = 80, L = 0) bandhead states of the ground 
g(K = 0) and fi'HK = 0) bands. Regular sequences at higher energy correspond 
to J3"y 2 (K = 0), /3"y 4 (K = 0) bands, etc. 



space portrait, shown for p = 0.2 in Fig. [2] is similar to the 
Henon-Heiles system (HH) Ir24ll with regularity at low energy 
[panels (bi)-(b2)] and marked onset of chaos at higher energies 
[panels (t^Hbs)]. The chaotic component of the dynamics in- 
creases with p and maximizes at the spinodal point p* = 0.546. 
The dynamics changes profoundly in the coexistence region, 
shown for p = 0.65, 0.741 in Fig. and £ = 0.05 in Fig. 
As the local deformed minimum develops, robustly regular dy- 
namics attached to it appears. The trajectories form a single 
island and remain regular at energies well above the barrier 
height Vb, clearly separated from the surrounding chaotic en- 
vironment. As £ increases, the spherical minimum becomes 
shallower, the HH-like dynamics diminishes and disappears at 
the anti-spinodal point £** = 0.354. Regular motion prevails 



for £ > where the section landscape changes from a single 
to several regular islands. The dynamics is sensitive to local 
degeneracies of normal-modes, as can be seen by comparing 
Fig. 3(di) for£= 1.094, corresponding to ep = e y , with Fig. 3(ei) 
for £=1.1. 

The quantum manifestations of the above rich classical dy- 
namics can be studied via Peres lattices {jc,-, 2?,} B25I1 . Here 
E, are the energies of eigenstates |z) of the Hamiltonian and 
x, = y/2(i\hd\i)/N. The lattices can distinguish regular from 
irregular states by means of ordered patterns and disordered 
meshes of points, respectively lEM l26ll . The particular choice 
of Xi can associate the states with a given region in phase space 
through the classical-quantum correspondence f3-x<^> Xj, ob- 
tained by comparing with the expectation value of hd in the 




Figure 4: (Color online). Peres lattices {*/,£j} for N = 50, L = 0,2,3,4 eigenstates of H' mt (p = p c ) = fl'J^ = f c ) (2), with hj = 1, /3 = 1-35 [panel (a)] and 
additional collective terms, Eq. |4}, with c,//!2 = 1, involving 0(3), 0(5) and 0(6) rotations [panels (b), (c), and (d)]. To enhance visibility, a small energy shift 
AEl = 0.005L(X + 1) is added. The classical potential shown, is the same in all cases. Notice in panels (a)-(b)-(c), the well-developed rotational bands (K = 0, 
L = 0, 2, 4) and (K = 2, L = 2, 3, 4) formed by the regular states in the deformed phase, which are distorted in panel (d). 



static condensate lUlll . The Peres lattices for L-0 eigenstates of 
Hi nt © with N = 80, are shown on the bottom rows of Figs. [2]|3] 
overlayed on the classical potentials V(x,y = 0). Forp = 0, the 
Hamiltonian (l2ab has U(5) dynamical symmetry with a solvable 
spectrum Ei=2Ji2\fi\N-\ +(l-fiQ)n c i]n c i. For large N and replac- 
ing Xi by p, the Peres lattice coincides with V ! (p = 0), a trend 
seen in Fig. 2(a). Whenever a deformed minimum occurs in the 
potential, the Peres lattices exhibit regular sequences of states, 
localized in the region of the deformed well and persisting to 
energies well above the barrier. They are related to the regular 
islands in the Poincare sections and are well separated from the 
remaining states, which form disordered (chaotic) meshes of 
points at high energy. The number of such sequences is larger 
when the potential well is deeper. The regular L = states 
form bandheads of rotational sequences L = 0, 2, 4, . . . (K — 
bands). Additional /T-bands with L-K,K+ l,K + 2, ... can 
also be identified. An example of such regular K = 0, 2 bands 
for H[ nt at the critical point, is shown in Fig. 4(a). In the nu- 
clear physics terminology, the lowest K — band refers to the 
ground band and excited /T-bands correspond to multiple p 1 and 
y vibrations about the deformed shape with angular momen- 
tum projection K along the symmetry axis. The states in each 
band share a common intrinsic structure as indicated by their 
nearly equal values of (fid), and a coherent decomposition of 
their wave functions in the rotor basis [27]. The occurrence of 
such a pure ordered band structure amidst a complicated en- 
vironment, indicates the relevance, for QPTs, of an adiabatic 
separation of modes |28] and possibly partial symmetries J29ll . 
for a subset of states. 

So far the discussion involved the intrinsic part of the Hamil- 



tonian (HJ. The collective part, which does not affect V(J3, y) 
can be transcribed in the form [ 19] 

H co i = c 3 [ C (3) -6h d ]+c 5 [ C (5) - 4« d ] 

+ c 6 [C m -5N], (4) 

where Cq denotes the quadratic Casimir operator of the group 
G, as defined injjjl] and c, = a/N(N - 1). The kinetic 
0(3), 0(5) and 0(6) terms involve collective rotations asso- 
ciated with the Euler angles, y and j3 degrees of freedom, re- 
spectively. Fig. 3] shows the Peres lattices corresponding to 
L — 0, 2, 3, 4 eigenstates of H m t at the critical-point, plus added 
rotational terms one at a time. As seen in Figs. 4(b)-4(c), the 
C3 and C5 terms preserve the, previously mentioned, ordered K- 
bands of H m t, Fig. 4(a). The calculated spectrum of these bands 
resembles a rigid rotor [L(L +1) splitting] for the 0(3) term 
and a rotor with centrifugal stretching for the 0(5) term [22]. 
In contrast, the regular band-structure is strongly disrupted by 
the 0(6) term [Fig. 4(d)]. The latter couples the deformed and 
spherical configurations [22] and mixes strongly the regular and 
irregular states. Only the 0(6) rotations involve the motion 
in the /3 variable lfl9ll . highlighting the importance, in QPTs, 
of the coupling of the order parameter fluctuations with soft 
modes 113011 . These results demonstrate the advantage of using 
the resolution of the Hamiltonian (Q]i in studies of QPTs, since 
a strong 0(6) term in the collective part can obscure the simple 
patterns of the dynamics disclosed by the intrinsic part. 

In summary, we have presented a comprehensive analysis of 
the dynamics across a generic first order QPT between stable 
spherical and deformed configurations in the IBM framework. 
The intrinsic part of the Hamiltonian determines the Landau po- 
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tential, and its classical analysis reveals a change in the system 
from an AO- and HH-type of dynamics on the spherical side, 
into a pronounced regular dynamics on the deformed side of 
the transition. The dynamics inside the coexistence region is 
robustly regular and confined to the deformed well, in marked 
separation from the chaotic behavior ascribed to the spherical 
well. The coexistence of regular and chaotic motion persists 
in a broad energy range throughout the coexistence region and 
is absent outside it. This simple pattern manifests itself also 
in the quantum analysis, disclosing regular rotational bands in 
the deformed region, which persist to energies well above the 
barrier and retain their identity amidst a complicated environ- 
ment. These ramifications of a divided phase space structure 
are observed at any Bq > 0, but are more pronounced for higher 
barriers (larger fio) ll lll27ll . Kinetic terms in the collective part 
of the Hamiltonian involving rotations in the orientation (Euler 
angles) and triaxiality (y) variables, preserve the ordered band- 
structure, while collective rotations in the deformation (J3) vari- 
able can disrupt it by mixing regular and irregular states. Our 
results are of interest not only for nuclei, but also to other inter- 
acting systems undergoing a first order QPT They demonstrate 
a clear connection between order, chaos and structural changes 
of coexisting phases in such systems. 
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